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ABSTRACT 

Context. The solar-like pulsator HD 49385 was observed with the CoRoT satellite during 137 days. The analysis of 
Oh' its oscillation spectrum yielded precise estimates of the mode frequencies over nine radial orders and pointed out some 

unusual characteristics: there exist some modes outside the identified ridges in the echelle diagram and the curvature 
^0 ■ of the I = 1 ridge significantly differs from that of the I = ridge. 

Aims. We first search for stellar models reproducing the peculiar features of the oscillation spectrum of HD 49385. 
Having shown that they can be accounted for only by a low-frequency 1=1 avoided crossing, we investigate the 
information which is brought by the mixed modes about the structure of the core of HD 49385. 

Methods. We propose a toy-model to study the case of avoided crossings with a strong coupling between the p-mode 
and g-mode cavities in order to establish the presence of mixed modes in the spectrum of HD 49385. We then show 
that traditional optimization techniques are ill-suited for stars with mixed modes in avoided crossing. We propose a 
, new approach to the computation of grids of models which we apply to HD 49385. 

Results. The detection of mixed modes leads us to establish the post-main-sequence status of HD 49385. The mixed 
mode frequencies suggest a strong coupling between the p-mode and g-mode cavities. As a result, we show that the 
amount of core overshooting in HD 49385 is either very low (0 < a ov < 0.05) or moderate (0.18 < a ov < 0.20). The 
mixing length parameter is found to be significantly lower than the solar one (qcgm = 0.55 ± 0.04 compared to the 
solar value a® = 0.64). Finally, we show that the revised solar abundances of Asplund give a better agreement than 
the classical ones of Grevesse & Noels. At each step, we investigate the origin and meaning of these seismic diagnostics 
in terms of the physical structure of the star. 

Conclusions. The subgiant HD 49385 is the first star for which a thorough modeling was led, trying to reproduce all the 
properties of an avoided crossing. It gave the opportunity to show that the study of the coupling between the cavities 
in these stars can provide valuable contribution to the questions of core overshooting, the efficiency of convection and 
' the abundances of heavy elements in stars. 

Key words. Stars: oscillations - Stars: interiors - Stars: evolution - Stars: individual: HD 49385 
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Fig. 1. Evolution of the eigenfrequencies of I = 1 modes 
as a function of the age for a 1.3 M@ model. The frequen- 
cies are normalized by the square root of the mean density 
^/GM/R 3 , where M and R are the stellar mass and radius. 



(|Di Mauro et all I2004L ICarrier et all I2005D and to set an 
upper limit for the overshooting, but the data were not 
precise enough to further constrain the internal structure of 
the star. The recent development of space missions CoRoT 
(jBaglin et al.ll2006t ) and Kepler (|Koch et al.ll2010l) , by pro- 
viding long quasi uninterrupted time series of high precision 
photometric data, has opened new opportunities for the de- 
tection of mixed modes. 

In this paper, we investigate the case of the GO-type star 
HD 49385, which was observed with the satellite CoRoT 
over a period of 137 days between October 2007 and March 
2008. The analysis of the time seri es has shown that the 
star exhibits solar-like oscillations (Deheuve ls et al.l f2010. 
further referred to as D10) and the authors could unam- 
biguously identify modes of degree £ — 0,1,2 over nine 
radial orders. Precise estimates of the mode properties (fre- 
quencies, linewidths, amplitudes) were obtained by fitting 
Lorentzian profiles on the observed spectrum. The authors 
also pointed out several striking features of the oscillation 
spectrum of HD 49385. They detected significant peaks 
which do not follow the expected pattern of high-radial- 
order p modes. Based on the low value of the observed sur- 
face gravity (log g = 4.00 ± 0.06, D10) which indicates that 
the star is certainly an evolved object, they mentioned the 
possibility that some of these peaks might be the signature 
of mixed modes. The analysis of the oscillation spectrum 
also showed that the curvature of the £ = 1 ridge in the 
echelle diagram unexpectedly differs from that of the £ — 
ridge at low frequency. 

We here perform a modeling of HD 49385 based on the 
spectroscopic and seismic constraints derived for the star 
by D10. We present in Sect. [2] a preliminary modeling of 
the star which shows that main-sequence models fail to re- 
produce the peculiar curvature of the observed £ = 1 ridge. 
We then investigate the possibility that this feature might 
be caused by the presence of mixed modes in the oscillation 
spectrum of HD 49385. Until now, theoretical studies led 
on avoided crossings all made the assumption that only two 
modes are involved a nd they neglected the contributions of 
the other modes (e.g. IChristensen-Dalsgaardl[l98ll . iGabriell 



Table 1. Fundamental parameters of HD 49385 measured 
from spectroscopic and photometric observations by D10 
(top) and estimated using the observed value of Av and 

I'max {bottom). 



Measured Parameters 


T cff (K) 


6095 ± 65 


(logp) spec tro 


4.00 ± 0.06 


[Z/X] (dex) w 


+0.09 ± 0.05 


log(L/L ) 


0.67 ±0.05 


Parameters estimated from Av and ^ max 


M/M & 


1.31 ±0.12 


R/Rq 


1.96 ±0.07 


(logp) sc ismo 


3.97 ±0.02 



Notes. The metallicity is defined as [Z/X] = 

\og[(Z/X)/(Z/X) e ]. 

Il980f ). We study in Sect. [3] the case of avoided crossings 
when the coupling between the p-mode cavity and the g- 
mode cavity is too strong to consider this phenomenon as a 
two-mode only interaction and we show that they can gen- 
erate a distortion of the ridge comparable to the observed 
one. This leads us to show that the seismic properties of 
HD 49385 can only be accounted for by the existence of an 
I = 1 avoided crossing. We then point out in Sect.|4]the lim- 
itations of traditional modeling techniques when studying 
stars with avoided crossings and we propose a method to 
remedy this. This method is applied to search for optimal 
models of HD 49385 in Sect. [5J with the aim to determine 
the constraints that the frequencies of the mixed modes 
bring on the structure of the inner regions of HD 49385. In 
Sect.|H we discuss these results and investigate their mean- 
ing in terms of internal structure and physical processes. 

2. Characteristics of HD 49385 and first step 
modeling 

2.1. Observational constraints 

We first give a brief overview of the observational con- 
straints which were derived for this star in previous studies. 

2.1.1. Surface constraints 

The surface observables of HD 49385 were derived by 
D10 based on a detailed analysis of two high-quality spec- 
tra obtained with the NARVAL spectrograph and on the 
Hipparcos measurement of the star's parallax. We here use 
their results which are recalled in Table [TJ 

2.1.2. Seismic constraints 

The star HD 49385 is a solar-like pulsator which has been 
observed with the CoRoT satellite during 137 days. The 
ridges of degrees ^ £ ^ 2 were unambiguously identified 
in the echelle diagram of the power spectrum (D10). The 
frequencies of these modes were determined over nine radial 
orders by fitting Lorentzian profiles on the detected modes, 
using a maximum likelihood estimation. 

The star has a mean large separation of (A^) b s = 
56.3 ± 0.5 fiRz and the frequency at the maximum of the 
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Fig. 2. Profile of the I — (squares and blue line) and 
£ = 1 (triangles and red line) observed large separation 
represented with 1-a error bars. The blue dotted line shows 
the result of the least-square-fit to a sinusoid of the I = 
large separation (see text). 



signal is v max = 1010 ± 10 /xHz (the error bars on these 
m easures were determined based on the method prescribed 
bylM 

osser fc A ppourchaux 2009J) . These values can be used 
to derive a first esti mate of the stellar mass and radius us- 
ing scaling relations. iGai et all (|201 lh recently noted that 
since the errors in these estimates of the mass and radius 
are positively correlated, the error we obtain on the quan- 
tity M/R 2 and thus on log g is small. We reported in Table 
[I] the estimated values of M, R and log g which will later 
be confronted to the ones of our best fitting models. We 
notice that the value we obtain for log<? is consistent with 
the spectroscopic measure of this quantity and has a much 
lower error, as expected. 

A clear oscillation can be detected in the profile of the 
£ = large separation, as can be seen in Fi g. [2] Following 
the procedure proposed bv I Roxburgh! (|2009t ). we estimated 
a period of 230 ± 30 /iHz for this oscillation. This indicates 
an acoustic depth of the glitch responsible for the oscillation 
of either r g i itch /r* = 0.76 ± 0.04 or r glitch /r^ = 0.24 ± 0.04, 
where t* is the total acoustic radius of the star (we recall 
that two glitches at acoustic depths of t and — r result in 
an oscillation of the eigenfrequencies with the same period) . 

The analysis of the oscillation spectrum of HD 49385 
led D10 to remark several unaccounted for features. First, 
contrary to what is expected for main sequence solar-like 
pulsators, the profile of the £ — 1 large separation signifi- 
cantly differs from that of the £ = large separation, es- 
pecially in the low-frequency part of the observations (see 
Fig. [2]). In an echelle diagram, this translates into a growing 
difference between the curvatures of the £ — and £ = 1 
ridges at low frequency. Secondly, several peaks were found 
to be significant even though they lie outside of the iden- 
tified ridges. One of them, the peak labeled as tt\ by the 
authors, has a posterior probability of being due to noise 
as low as 10 -5 . It lies close to the £ = ridge but not in its 
direct continuity and it would cause an abrupt step in the 
£ = large separation profile if it were identified as a radial 
mode. One of the challenges of the modeling of HD 49385 
will be to understand these peculiarities. We address this 
question in the following section. 



2.2. Preliminary modeling of HD 49385 

2.2.1. Properties of the models 

All the mo dels are comp uted with the stellar evolution code 
CESAM2K (jMorell [19971 ) and the mode frequencies are de- 
rived from them usi ng the Liege Oscillation Code (LOSC, 
iScuflaire et al.|[2008h . We used the OPAL equation of stat e 
and opacity tables as described in iLebreton et al.l (|2008l) . 
The nuclear r eaction rates are com puted using the NACRE 
compilation ([Angulo et al.l [l999h . The atmosphere is de- 
scribed by Eddington's grey law and is connected to the 
envelope at an optical depth of r = 10 to ensure the v alid- 
ity of the diffusion approximation (|Morel et al.l [1994). We 
computed models using alternately the mixture of heavy 
elements of iGrevesse fc Noelsl (| 19931) , further referr e d to a s 
GN93, and the revised mixture of lAsplund et al.l (|2005T ) . 
referred to as AGS05. 

The convective regions are treated usin g the Canuto- 
Gold man-Mazzitelli (CGM) formalism (Ca nuto et al.l 
1996). It involves a free parameter, the mixing length, 
described as a fraction «cgm of the pressur e scale height. 
A cal ibration on the Sun gives «cgm = 0.64 (jSamadi et al.l 
120061 ). In this work, the mixing length will be considered 
as a free parameter. 

The radius R cc of the convective core is determined by 
the Schwarzschild criterion. Overshooting can be included 
in the models as an extension of the motion of convective 
eddies over a distance d ov outside the core, which is defined 
as 

d ov = a ov min (i? cc , H P ) (1) 

where Hp is the pressure scale height. In this study, we 
consider an instantaneous mixing in the overshooting re- 
gion, meaning we assume that the time scale of the mixing 
is much smaller than the time scale of the evolution of the 
star. In this case the overshooting zone is fully mixed. The 
temperature gradient is taken to be equal to the adiabatic 
gradient. 

Our models are computed neglecting microscopic diffu- 
sion, in order to limit the computational time of our grids 
of models. However, the effect of microscopic diffusion on 
our conclusions is studied in Sect. [51 

2.2.2. Comparison criterion between models and observations 

As is usually done, we compare the stellar models and the 
observations by computing the merit function \ 2 defined as 

r /nmod /nobs 1 ^ 

X 2 -E 1 ~^ ( 2 ) 
i=i 1 4 J 

where 0° hs , i = 1, N represent the N observables selected 
to constrain the models, cr° bs their error bars and Of 10 ' 1 the 
values of these parameters for the computed models. 

To constrain the models, we used the classical observ- 
ables T c ff and L/Lq. We note that the measured value 
of log<7 could also be used. However its value is not well 
enough constrained by spectroscopic measurements to al- 
low it to discriminate among the different models and we 
will therefore no longer mention it in the following. The set 
of observables is completed by the estimates of the p-mode 
frequencies of degrees ^ £ ^2 obtained by D10. We note 
that the authors also mentioned the possible detection of 



4 



S. Deheuvels and E. Michel: Constraints on the structure of the core of HD 49385 via mixed modes 



two £ = 3 modes (modes ~ki and 7^) and proposed estimates 
of their frequencies. However, given their weak signal-to- 
noise ratio, we chose not to include them in our set of ob- 
servables. Additionally, one of the I = 2 modes was found to 
overlap the closest radial mode (vo,ii = 855.3 /iHz). Since 
none of all the models we computed in this study present 
this particularity, we supposed that this phenomenon is due 
to the lower signal-to-noise ratio of £ = 2 modes around 
the edges of the frequency domain of the oscillations. This 
1 = 1 mode was therefore not included among the observ- 
ables <D° hs . 

To compare the eigenfrequencies between models 
and observations, we know that it is necessary to 
correct them from the effects of our improper mod- 
eling of the structure in the surfa c e layers (e.g. 
IChristensen-Dalsgaard fc Thompsonlll997r l. iKieldsen et al.l 
(2008) found that in the case of the Sun, the differences be- 
tween the observed frequencies and those of the best solar 
model were well approximated by a power law of the form 

/ obs\ b 

<|-<T = af^l (3) 

where i>q is the frequency of the maximum of the signal. 
The authors tested this law on a few other solar-like pul- 
sators and found a reasonable agreement. For this purpose, 
they assumed that the exponent found in the solar case 
(b = 4.90) could also be used for other targets. We ap- 
plied a correction of this type on the frequencies of our 
models. Only our cas e is somewhat different from that of 
IKieldsen et all ([20081) : our models are computed using the 
C ESAM2K code and we treat con vection with the forma lism 
of ICanuto etTai] (| 19961) whereas IKieldsen et al.l 1)20081 ) use 
ASTEC and describe convection using the traditional mix- 
ing length theory. It would therefore be irrelevant to use the 
same exponent as the one found by the authors in the solar 
case. We computed a sola r model mat c hing t he set of solar 
eigenfrequencies found bv iGellv et a l. ( 2002). We then fit- 
ted the power law given by Eq.|3]on the differences between 
the observed frequencies and those of our solar model. We 
obtained a somewhat smaller exponent (6 = 4.25) which 
is used to correct the frequencies in the following. We also 
note that the modes which have a mixed behavior are ex- 
pected to be less sensitive to the surface effects. To take this 
into account, for non-radial modes, the surface correction 
given by Eq. [3] was multiplied by a factor Q~\, where Q n ,l 
corresponds to the ratio between the inertia of the mode 
and the inertia of the closest radial mode, as prescribed by 
lAerts et al.l (|2010h . 
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Fig. 3. Evolutionary tracks of two models fitting the posi- 
tion of HD 49385 in the HR diagram. The observed values 
of the effective temperature and luminosity of HD 49385 
are indicated within 1-a error bars by the box. The blue 
dashed line corresponds to a MS model, and the red solid 
line to a PoMS model. 
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Fig. 4. Echelle diagram of the power spectrum of HD 49385 
derived from 137 days of CoRoT data, folded with a mean 
large separation Av = 56.3 /iHz. The frequencies of a MS 
model and a PoMS model are overplotted in blue and red, 
respectively. The models have been computed to reproduce 
the observed £ — large separation as well as the position 
of the star in the HR diagram and the PoMS has an £ — 1 
avoided crossing at low frequency. Squares represent i = 
modes, triangles £ — 1 modes and diamonds £ = 2 modes. 



2.2.3. Evolutionary status of HD 49385 



The low value of the spectroscopic log<? obtained for 
HD 49385 suggests that the star is evolved. We can find 
models both in the main sequence (hereafter MS) and post 
main sequence (PoMS) stage, which fit the position of 
HD 49385 in the Hertzsprung-Russell (HR) diagram (see 
Fig. [3]) . We are facing a problem which is quite common- 
place for evolved objects: a degeneracy of MS models and 
PoMS models which results in an uncertainty on the evo- 
lution ary stage of the studied object (see e. g. Procyon 



A in iBarban et al.l Il999l iProvost et all I2006L rj Boo in 



iDi Mauro et al.l l20031 



The internal structure of these two families of models is 
very different. The main sequence models are close to the 
exhaustion of their hydrogen reserves in the center (X c ~ 
0.1) and they have a small convective core spreading over 
about 5% of the stellar radius. The post main sequence 
models are burning hydrogen in a thin layer located above 
the limit of the convective core which existed during the 
main sequence stage. There are no longer any convective 
regions in the center of these models. 
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2.2.4. Main sequence models 

We computed a grid of main sequence models with varying 
masses, ages, helium abundances, mixing length parameters 
and metallicities. We found models which fit quite well the 
observed variations of the £ = large separation and the 
position of the star in the HR diagram. The echelle diagram 
of one of them is represented in Fig. 2) However, none of 
the computed MS models are able to reproduce the peculiar 
curvature of the I = 1 ridge which we mentioned before. 
The values of the x 2 function we obtain for MS models 
are all above 2000. This yields a reduced \ 2 above 100, 
which indicates a very poor match with the observations. 
The high value of \ 2 f° r MS models is almost entirely due 
to the contribution of the I = 1 modes at low frequency 
(the four £ = 1 modes with the lowest frequencies account 
for more than 90% of the x 2 value for the best MS models). 

Besides, no MS model can explain the presence of the 
peak 7Ti. It cannot correspond to the signature of a mixed 
mode, since MS models which reproduce the surface observ- 
ables of HD 49385 are not evolved enough to have mixed 
modes in the frequency domain of the observations. This 
peak could not either be identified as an £ = mode since 
the closest radial mode in the models lies about 5 fiKz away 
from it (i.e. at more than 20a). 



2.2.5. Post main sequence models 

Contrary to MS models, some PoMS models which fit the 
surface parameters of HD 49385 have £ = 1 mixed modes 
in the frequency domain of the observations. We have con- 
sidered the possibility that the difference we observe be- 
tween the curvatures of the I = 1 ridge and the £ = 
ridge might be caused by an I — 1 avoided crossing. So far, 
the studies which were led about avoided crossings made 
the assumption that only two modes were involved and 
they neglected the contributions of the o ther modes (e.g. 
IChristensen-Dalsgaardll98"li [Ga briel 1983). In the next sec- 
tion, we investigate the case when the coupling between the 
cavities is such that this hypothesis is no longer valid and 
we try to determine the effect it has on the curvature of the 
ridge. 




Fig. 5. Variations of the eigenfrequencies ui± of the system 
with parameter A, for different values of the coupling term 
a. 



sponding to the following system of equations: 



dt 2 

dt 2 



-wi(A) 2 yi + ay 2 
-w 2 (A) 2 y 2 + ayi 



(4) 



where a is the coupling term between the two oscillators 
and wi(A), u>2(A) are the eigenfrequencies of the uncoupled 
oscillators (particular case where a = 0). They are chosen 
to depend on a parameter A, used to model the change 
of the dimensions of the cavities as the star evolves. We 
suppose that for a certain A = Ao, the frequencies of the 
uncoupled oscillators cross, i.e. wi(Ao) = cj2(Ao) = loq. 
Solving Eq. 2J we obtain the two following solutions: 



± 



1 



,2^2 



4a 2 



(5) 



If the coupling term a is very small compared to the differ- 



ence between the eigenfrequencies (a <C 



\w\- 



;|), then the 



eigenfrequencies of the system are close to u>i and u>2- If, on 



the contrary, 
can be approximated by 



<C a i,2, then the eigenfrequencies 



3. Avoided crossings with strong coupling 

3.1. Analogy with harmonic oscillators 

We here develo p an analogy introduced by 
iDeheuvels fc Michell (|2010h to gain insights on the 
characteristics of avoided crossings involving n coupled 
modes. We recall its main points, for the sake of clarity, 
and we add some details about the choice of the coupling 
term. 



«4 



UJn ± a. 



(6) 



To analyze these solutions, we choose the eigenfrequencies 
of the uncoupled oscillators a>i(A) and w 2 (A) such that they 
simulate a p mode and a g mode. By reference to Fig. [Tl we 
impose uji(X) = 1 and W2(A) = A. In that case, an avoided 
crossing occurs in the system around Ao = 1. The varia- 
tions of the two eigenfrequencies w±(A) are shown in Fig.[5J 
The two modes clearly exchange natures during the avoided 
crossing and the intensity of the phenomenon depends on 
the strength of the coupling between the oscillators. 



3.1.1. Avoided crossing with two modes 

Christensen-Dalsgaardl (120031) . based on the work of 
von Neuman fc Wignerl (|l929ft . proposed a simple analogy 
which captures the main aspects of avoided crossings be- 
tween two modes. He considers the two cavities of the star 
as two coupled harmonic oscillators y\{t) and y 2 (£)i re ~ 



3.1.2. Avoided crossing with n modes 

The previous analogy relies on the fact that two modes only 
arc affected during an avoided crossing: this corresponds to 
neglecting the coupling between the two considered modes 
and the other modes in the spectrum. Let us now push the 
analogy a step further, and consider the case where these 
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Fig. 6. Variations of the eigenfrequencies of (n— 1) p modes 
coupled to a j mode which undergoes avoided crossings 
with the p modes (here, n = 10). The dashed lines corre- 
spond to a "weak coupling" (7 = 0.01), and the full lines 
to a "strong coupling" (7 = 0.06). 



these coupling terms in our analogy. As A increases, the 
g mode will experience successive avoided crossings with 
all the p modes. According to Eq. [51 the deviation of the 
eigenfrequency very close to the i th avoided crossing is 
Slo ~ ai/(2coi). To ensure that all the avoided crossings 
have a comparable intensity (as seems to be roughly the 
case in the models, see Fig.QJ, we choose the different cou- 
pling terms ai such that 

oti = 7Wj (9) 

where 7 describes the strength of the coupling between the 
p-mode cavity and the g-mode cavity. We obtain the fol- 
lowing system of equations 



dV(t) 
dt 2 



d 2 y n -i{t) 
dt 2 

d 2 yn(t) 

dt 2 



(10) 
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Fig. 7. Echelle diagrams of the eigenfrequencies of 10 cou- 
pled harmonic oscillators, at a given "time" A in the vicinity 
of an avoided crossing. The left panel presents the case of 
a weak coupling (7 = 0.01) and the right panel the case 
of a strong coupling (7 = 0.06). The frequencies which are 
circled correspond to the ones which are used to estimate 
the strength of the coupling in Sect. |6j 

other coupling terms play a significant role. We consider n 
oscillators, instead of two only. We choose n — 1 of them to 
simulate high-radial-order p modes, by giving them equidis- 
tant eigenfrequencies which are constant with A. The last 
oscillator simulates a g mode. We have taken 



uJi(A) = Wj_i(A) + Aw, 
w„(A) = A 



for 



1 



(7) 
(8) 



The numerical values of wi(A) and Aw (1 and 0.1 respec- 
tively) have been chosen to roughly match Fig. [1] 

We then introduce a coupling ai between the g mode 
and the i th p mode. Normally the p modes should also 
be coupled to each other. However, since their eigenfre- 
quencies remain equidistant at all times, we can neglect 



By writing the different oscillators Ui(t) — c; exp(— iuit), 
the eigenfrequencies of the system are found by solving the 
eigenvalue problem AC = ui 2 C with 



A 



w? 

















-7W1 



W 



n-l 
7^n- 



-7W1 
-7W2 

-7 w n-l 



(11) 



and C = [ci, • • • ,c n ]. 

The solution is plotted in Fig. [6] for two different val- 
ues of 7. For each value of the parameter A, we can plot 
an echelle diagram with the eigenfrequencies of the coupled 
system (Fig. [7]). In the case of a weak coupling, the ridge re- 
mains straight and the approximation that only two modes 
are affected by the avoided crossing is legitimate. However, 
for a strong coupling, more p modes in the neighborhood of 
the avoided crossing have a mixed behavior and the ridge 
is distorted as shown in Fig. [JJ 

3.2. Comparison with stellar models 

In a star, the coupling is due to the existence of an evanes- 
cent zone between the two cavities. The wider this zone, the 
weaker the coupling. Since the Lamb frequency increases 
with the dergee £, the coupling is stronger for modes of 
small degrees. We represented in Fig. [8^i the propagation 
diagram of a 1.3 M© model, evolved enough so that avoided 
crossings occur for low-degree modes (post main sequence) . 
For this model, the ridges of the modes of degrees 1 ^ £ ^ 3 
are shown in an echelle diagram in Fig. [S}d. Clearly, for 
modes of degree t ^ 2, the currently adopted hypothesis 
that only two modes are affected by the avoided crossing is 
valid, the rest of the ridge being almost unaltered. However, 
for the £ — 1 avoided crossing, the coupling is strong enough 
so we observe a distortion in the ridge which is very compa- 
rable to the one we obtained in our simple analogy (bottom 
panel of Fig. EJ). 
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Fig. 8. Left: Propagation diagram of a 1.3 M© post main sequence model. The Brunt- Valsala frequency is represented by 
the full black line and the Lamb frequencies by the dashed lines {1 = 1 red, I = 2 green, £ = 3 purple). The propagation 
regions of three mixed modes of degrees £ = 1,2 and 3 are represented by the full colored lines and their evanescent 
regions by the dotted colored lines. Right: Echelle diagram of the same model {£ = blue squares, 1=1 red triangles, 
1 = 2 green diamonds, i = 3 purple crosses). The mixed modes represented in the left panel are circled. 



3.3. Application to the case of HD 49385 

By fine-tuning the mass and age of PoMS models, we man- 
aged to find a model with an £ = 1 avoided crossing in 
the low-frequency part of the frequency range of the ob- 
servations (around 750 fiRz) and which also reproduces the 
observed mean value of the large separation (we will present 
in Sect. |4] a method to find these models in a systematic 
way). The echelle diagram of this model is represented in 
Fig. |U We observe that the curvature of the 1=1 ridge 
is far better reproduced by this model than by MS mod- 
els. Interestingly, we note that by adjusting the age of the 
model, the frequency of the mode which behaves mainly as 
a g mode in the avoided crossing matches the frequency of 
the peak 7ti we mentioned before. 

We conclude that the observed oscillation spectrum of 
HD 49385 can be satisfactorily reproduced only if an £ = 1 
avoided crossing causes a distortion the 1 = 1 ridge in the 
echelle diagram. We thus obtain a firm detection of mixed 
modes in the spectrum of the star. We note that the avoided 
crossing occurs at the far bottom of the frequency range of 
the observations, which explains why we only detect the 
upper half of the ridge distortion described in Sect. El We 
have also been able to identify the mode which has a mainly 
g-mode behavior in the avoided crossing (corresponding to 
the peak tt±). This will be very valuable when conducting a 
thorough modeling of the star, since this mode is the most 
sensitive to the center of the star. Finally, the detection of 
an I = 1 avoided crossing in the spectrum of HD 49385 
enables us to establish that the star is in a PoMS stage. 



4. Stellar modeling using the avoided crossing: 
potential and strategy 

4.1. Potential of avoided crossings with strong coupling 

Mixed modes have a high potential in terms of seismic di- 
agnostics since they are sensitive to the structure of the 
most central regions of the stars, which remains poorly un- 
derstood. 



The frequencies of mixed modes depend on the profile of 
the Brunt- Valsala frequency in the stellar core. To under- 
stand the information which we can expect to derive from 
them, it is convenient to write the Brunt- Valsala frequency 
as a function of the following gradients: 



<91nT 



ad 



ainP 



91nT 
dlnP 



V/i 



d In fi 
dlnP' 



(12) 



where P, T and /i correspond to the stellar pressure, tem- 
perature and mean molecular weight, respectively. The sub- 
script S indicates that the definition is valid for constant 
entropy The Brunt- Valsala frequency can be split in two 
different contributions: 



N z = 



JJ_ 

Hp 



(V ad - V) 



H, 



-Vp, 



(13) 



where g is the gravity inside the star and Hp is the local 
pressure scale height. The first term in the right hand side 
of Eq. [T3] depends on the temperature stratification. The 
right term characterizes the dependance on the gradient 
of mean molecular weight V//. These two components are 
shown in Fig. [HI for the model presented in Sect. 13.21 

Two different features of avoided crossings with strong 
coupling bring information about the structure in the center 
of the star: 

1 . The frequency at which the avoided crossing occurs 
The knowledge of the frequency of an avoided crossing 
provides an estimate of the frequency of the g mode 
which is involved. This opportunity is extremely pre- 
cious, since the observation of pure g modes in main 
sequence stars is currently out of reach. In t he frame- 
work of the asymptotic theory (iTassoull 198.3) , the fre- 
quency of a g mode of radial order n and degree £ can 
be approximated by 



(n — l/2j7r 



(14) 



where r\ and ri are the mode turning points in 
the g-mode cavity. In the following, we will consider 
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Fig. 9. Propagation diagram for a 1.3 Mq post main se- 
quence model which presents an I = 1 avoided crossing. 
The Brunt- Vaisaia frequency (black solid line) is split in 
two contributions: the part linked to the temperature strat- 
ification (long-dashed blue line) and the one linked to the 
gradient of mean molecular weight (red solid line). The 
Lamb frequency for £ — 1 modes is represented by the 
black dashed line. The shaded area indicates the evanes- 
cent zone. The radii r± and T2 correspond to the turning 
points of the g-mode cavity and r% the turning point of the 
p-mode cavity. 
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Fig. 10. Top: Evolution of the frequencies of I = modes 
(dashed black lines) and I — 1 modes (full black lines) with 
age for a 1.3 Mq model. The frequencies of the first two 
I — 1 uncoupled g modes are also represented (gi mode in 
red and g2 mode in green) . Bottom: Evolution of the central 
density p c with age for the same model. 



avoided crossings involving low-radial-order g modes, 
for which the asymptotic approximation is no longer 
valid. However, their frequencies mainly depend on 
/J" 2 Ndr/r. Having access to the frequency of a g mode 
gives information about the structure inside the g-mode 
cavity (r\ ^ r ^ T2 in Fig. [§]). In particular, it is 
sensitive to the peak left in the profile of V/x by the 
withdrawal of the convective core, for stars massive 
enough to have had one during the main sequence. 
This is the reason why the detection of mixed modes 
is expected as a means of constraining the amount of 
overshooting at the boundary of c onvective cores (e.g. 
iDziembowski Sz Pamvatnvkhlll991f) . 

2. The intensity of ridge distortion 

We showed in Sect. [3] that more than two modes are 
affected by I = 1 avoided crossings. The modification 
of their frequencies (characterized by a distortion of 
the £ = 1 ridge in the echelle diagram) depends on 
the strength of the coupling between the p-mode cav- 
ity and the g-mode cavity. This coupling mainly de- 
pends on the profile of the Brunt- Vaisaia frequency in 
the evanescent zone: the higher N(r) in this region, the 
stronger the coupling. The intensity of the ridge distor- 
tion should therefore bring information about the struc- 
ture of the star in the evanescent zone (i.e. in the region 
ri ^ T ^ ?"3 in Fig. 0, complementary to those given 
by the frequency of the avoided crossing. We will illus- 
trate this complementarity with the modeling of the star 
HD 49385 in Sect. U and M 



4.2. Limitations of traditional optimization techniques to 
model stars with avoided crossings 

As mentioned in Sect. [TJ numerous theoretical works have 
been led about avoided crossings. However, very few studies 
have tried to find stellar models fitting the properties of an 
avoided crossing. The first reason for this is that until very 
recently, the spectra in which the detection of mixed modes 
was claimed were too imprecise to lead such an investiga- 
tion. The second reason is that it presents some inherent dif- 
ficulties linked to the fact that the avoided crossings occur 
on a very short timescale (typically of the order of 1 Myr or 
less) compared to the stellar evolution timescale. The usual 
approach to model a star consists in computing a grid of 
models with varying stellar parameters (mass, age, etc.) to 
find the optimal model which minimizes the \ 2 function 
defined in Eq. [2l Applying this procedure with a time step 
of the order of the avoided crossing timescale is not feasible, 
since it would require to compute a tremendous number of 
models. On the other hand, with a larger time step, we have 
a very low probability to find models which correctly repro- 
duce the frequency of the avoided crossing and we therefore 
miss the best-matching models. We note that it would also 
be impractical to bluntly apply an automatic minimization 
using for instanc e the algorithm of Levenber g-Marquardt, 
as prescribed bv iMiglio fc Montalba n (2005). Indeed, this 
technique is based on the computation of the Hessian ma- 
trix at each iteration and it thus requires to input a step for 
the time derivatives. The gap between the stellar evolution 
timescale and that of the avoided crossing makes it impos- 
sible to find a satisfactory time step for the procedure. 
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4.3. Narrowing down the dimension of the model space 

We now try to adapt the traditional grid-of-model approach 
to the special case of targets for which mixed modes are 
detected. For this purpose, we show that we can overcome 
the obstacles we just mentioned by no longer considering 
the stellar age and mass as free parameters of our fit. 

4.3.1. Relation between f crc ,ss and the stellar age 

The frequency at which an avoided crossing occurs (further 
noted Across) corresponds to the frequency of the uncoupled 
g mode it involves and is therefore linked to the profile of 
the Brunt- Vaisala frequency in the core (see Sect. 14. Fig. 
[TU1 shows the variations of the frequencies of the first two 
uncoupled g modes during the evolution of a 1.3 M© model 
(the sharp variations of the mode eigenfrequencies around 
the age of 4050 Gyr in this plot correspond to the transition 
between the MS stage and the PoMS stage). 

During the main sequence, the first term in the right 
hand side of Eq. [T3] nearly vanishes in the center since 
the temperature gradient is almost adiabatic in the convec- 
tive core. The g-mode cavity essentially corresponds to the 
chemically inhomogeneous zone left by the withdrawal of 
the core. At this point, the Brunt- Vaisala frequency is still 
low, so that low-degree g modes do not reach the frequency 
range of high-order p modes and no avoided crossing can 
be observed. 

At the end of the main sequence, when the hydrogen 
reserves are exhausted in the core, the star is left with an 
almost isothermal helium core. The inner regions contract 
until nuclear reactions are triggered in the layers where hy- 
drogen remains. The evolution of the Brunt- Vaisala fre- 
quency is then mainly determined by the value of the cen- 
tral density p c . To explain this, we remark from Ea. ll3l that 
N depends on the factor g/ Hp, which can also be written 
pg 2 1 P. In the most central regions, the gravity term can be 
approximated as g(r) ~ Gp c r. In that case, we have 

■f ~ P ~ fV (15) 

tip f c lc 

where P c and T c are the central pressure and tempera- 
ture. The models show that the central temperature hardly 
changes in the isothermal core during the evolution in the 
post main sequence phase. The mean molecular weight is 
maintained constant since the chemical composition in the 
core does not vary. On the contrary, the contraction of the 
star causes the central density to increase by more than one 
order of magnitude compared to its value at the end of the 
MS, as can be seen in the lowest panel of Fig. [TU] As a re- 
sult, the factor g/Hp and thus the Brunt- Vaisala frequency 
mainly depend on the evolution of the central density p c . 
This is confirmed by Fig. ITOl which shows the tight relation 
between the evolution of Across and p c for a 1.3 Mq model. 

We conclude that the avoided crossing frequencies 
monotonically increase during the evolution of the star. We 
can therefore find for every model (with fixed mass and 
physics) one stellar age such that the avoided crossing oc- 
curs at the same frequency as in the observations. For the 
model represented in Fig. [TU1 an age of 4374 ± 0.8 Myr 
satisfies this condition (the error bar is determined by re- 
quiring a fit of v°^o SS within a 1-ct error bar) . The precision 
we obtain is such that we can consider the age to be en- 
tirely determined by the value of Across (for models at given 
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Fig. 12. Frequency of the avoided crossing for the models of 
Fig. llll which reproduce the observed mean large separation 
of HD 49385. The dashed line indicates the frequency of the 
observed avoided crossing v°^ ss for HD 49385. 

mass and physics) and we can cease to consider the age as 
a free parameter of the fit. This is very interesting since it 
solves the problem of the choice of a time step for the grid 
of models. 

We note that by using this method, we overlook the 
models which do not fit the frequency of the avoided cross- 
ing within 1-ct error bars. However, since avoided crossings 
occur on a small timescale compared to the stellar evolution 
timescale, modifying is cross in a model will result in almost 
no change in the other stellar parameters. It will only in- 
duce an increase in the contribution of the term including 
Across i n the x 2 defined in Eq.[5]with no significant change in 
the other contributions. Overlooking these models should 
therefore not modify the optimal model we obtain nor the 
error bars we derive from it. 

4.3.2. Correspondence between (Ai/, z/ cross ) and (M, age) 

We just saw that for a model with given mass and physics, 
we can find an age such that Across = ^°ross- However, it is 
important to notice that the mean large separation of the 
star monotonically decreases with the age since the stellar 
radius keeps increasing (even during the sharp transition 
between the main sequence and the post main sequence 
stage). For every model, there exists an age at which the 
observed value of the mean large separation is reproduced. 
This age is a priori different from the one for which the 
frequency of the avoided crossing is reproduced. So which 
are the models that fit both conditions simultaneously? We 
will further refer to this joined condition as condition C, i.e. 

_ J(A^) mod = (A^) bs ri(n 

U • 1 mod _ obs \ 1U ) 

L cross cross 

To answer this question, we search among the models 
fitting the observed value of the mean large separation those 
which also reproduce the frequency of the avoided crossing. 
As explained in Appendix [AJ the iso- Az/ region in the HR 
diagram roughly corresponds to a line of slope 5 (L oc T c 5 ff ). 
Fig. llll shows the evolutionary tracks of several models with 
the same given physics and different masses. We observe 
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Fig. 11. Left: Evolutionary tracks in the HR diagram for models of masses ranging from 1.23 M© (lightest symbols) to 
1.31 Mq (darkest symbols), with given physics. The filled circles indicate the points where the models reach (Av) = 
(A^) bs for HD 49385 (this location is materialized by the dotted line). Right: Evolution of the central density for the 
same models. The value of the mean large separation has been used as an indicator of the age, for more convenience. 



that as the stellar mass increases, the line of the Terminal 
Age Main Sequence (TAMS) gets closer to the line of iso- 
Av. Consequently, the more massive stars are closer to the 
TAMS when they reach the observed value of the large sep- 
aration. This means that they have smaller central density, 
as can be seen in Fig. [IT] (where we found it convenient 
to use the mean large separation instead of the age as an 
indicator of the evolution). Since we mentioned that the 
evolution of Across is determined by that of p c , we conclude 
that Across should decrease with increasing mass. This is 
confirmed by Fig. 1121 This result should be kept in mind, 
as it will frequently be used in the following. 

We reach the conclusion that there exists one and only 
one value of the stellar mass and age for which both the 
mean value of the large separation and the frequency of 
the avoided crossing are correctly fitted (for instance with 
the physics used in Fig. [11] and [T^J it corresponds to a 
mass of about 1.28 Mq). We will further note M and f 
the stellar mass and age which satisfy condition C. It is 
striking to remark that for a given physics, the two seismic 
parameters (Av) and v CTOSS alone provide an estimate of the 
stellar mass and age, without having to resort to the use of 
any classical constraints or other seismic parameters. This 
estimate is also extremely precise. For example, if we fix all 
stellar parameters other than the mass and age in the case 
shown in Fig. 1121 the uncertainty we obtain on the mass 
is as low as 4 x 10~ 4 M©. We therefore propose to narrow 
down the dimension of the model space by eliminating the 
mass and age from the set of free parameters, these two 
quantities being determined by ((Av) , v CIOSS ) . Of course, 
the values of M and f depend on the physics we use to 
model the star. Their dependance with the different stellar 
parameters will be studied in Sect. 14.51 after describing the 
method we adopted in this study to find the models which 
satisfy condition C. 

4.4. Searching for models satisfying condition C 

We now briefly describe the procedure we followed to obtain 
the values of M and f for each set of free parameters. 



As mentioned in Sect. 14.3.11 we consider the stellar age 
to be fixed by imposing v cross = Across- Until now, we have 
assumed that the frequency of the avoided crossing corre- 
sponds to the frequency of the uncoupled g mode it involves. 
Of course, when computing the eigenfrequencies of a model, 
we do not have access to this quantity. We thus chose to 
assimilate the frequency of the avoided crossing to the fre- 
quency of the mode which behaves mainly as a g mode. 
In the oscillation spectrum of HD 49385, it corresponds to 
the peak labeled as m, whose frequency has been found 
to be tVi = 748.6 ± 0.23 //Hz by D10, and in the models 
it is the mode ^1,11. For each computed model, the age is 
determined so that we have fi,n = v 7ri . Special care has to 
be taken in this procedure since the evolution of the £ = 1 
mode frequencies with the age is not monotonic due to the 
avoided crossings, as can be clearly seen in Fig. [TU] As a 
result, there are several ages (three) for which (v\ t \\ = v vi ) 
and only one among them corresponds to the one we search 
for. 

We need to determine the stellar mass for which the 
models fitting the frequency of the avoided crossing also 
reproduce the observed large separation. This mass is the 
one that minimizes the merit function 

2 _ ((Ai/) mod - (A^) obs ) 2 

X{Av) = ~2 \ LI ) 

where (Ai/) b s is the mean large separation of the observed 
radial modes, criAv) the 1-a error-bar and (Av) mo a the cor- 
responding quantity for the models (after correcting for 
the surface effects as described in Sect. I2.2.2[) . Only the 
radial modes are considered here since the large separa- 
tions of non-radial modes are affected by avoided crossings. 
We remark that for models reproducing ^°ross> the mean 
large separation of radial modes varies almost linearly with 
the stellar mass. The function X^au) * s therefore almost 
quadratic and the method of Newton is particularly effi- 
cient to determine the optimal mass M (obtained in only a 
couple of iterations). 
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Fig. 13. Variations of the mass M when changing the mixing length parameter (left), the initial helium abundance 
(middle) and the metallicity (right) from model So (whose properties are given in Table [2]). 
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Table 2. Parameters of the models used to determine the 
variations of the mass M (which verifies condition C) with 
the stellar parameters. 



Model 


Qconv 


Y 


Mixture 




Mass 


So 


0.56 


0.24 


AGS05 


Mo 


(1.263 M ) 


Sla 
Sib 


0.60 
0.60 






Mi 


M 
(1.295 M Q ) 


S2a 
S2b 




0.28 
0.28 




M 2 


Mo 
(1.202 M Q ) 


S3a 
S3b 






GN93 
GN93 


M 3 


Mo 
(1.308 M Q ) 



Fig. 14. Evolutionary tracks in the HR diagram of models 
S , Si a , Sib, S 2a , S2b, S 3a and S 3b (see text and TableHJ. 



4.5. Dependance of M on the stellar parameters 

By applying the procedure we just described, the mass M 
and age f can be determined for any given set of free pa- 
rameters. This gives us the opportunity to study the depen- 
dance of M with the different stellar parameters which are 
varied in this study (the mixing length parameter a conv , 
the initial helium abundance Yq, the metallicity [Z/X] and 
the amount of core overshooting a ov ). This will provide us 
with helpful insights when discussing the results of our op- 
timization in Sect. \E[ This should also give some hints of 
how the mass should be varied when studying other stars 
with avoided crossings. 

We start from a model Sq with a given set of stellar 
parameters (see Tabled]). From what has just been said, 
there exists one stellar mass Mq and age fo which verifiy 
condition C with these parameters. By varying the stellar 
parameters one by one we try to understand the variations 
of the optimal mass M with each of them. 

4.5.1. Dependance of M on the mixing length parameter 

We modified the mixing length parameter of model Sq with- 
out changing the other parameters and we determined the 



new mass M by applying the procedure described in Sect. 
14.41 Fig. IT3"k shows that M linearly increases with a conv . 

The increase of M with a conv can be understood as fol- 
lows. Let us modify the mixing length parameter from the 
value ao of model So to a value a\ > «o- If we keep the same 
mass Mo (model Si a in Tabled, an increase of the mixing 
length induces an increase of the convective flux in the enve- 
lope and therefore causes the star to contract slightly. Since 
the stellar luminosity remains unchanged, the decrease of 
the radius is compensated by an increase of the effective 
temperature. The evolutionary track is thus horizontally 
translated to the left in the HR diagram when the mix- 
ing length goes from ao to ot\ at constant mass (see Fig. 
113b ). As a result, for a conv = ot\, the star is further away 
from the TAMS when it reaches the observed value of the 
large separation and the frequency of the avoided crossing 
is too high compared to the observations. Based on what 
was said in Sect. I4.3T21 we need to increase the stellar mass 
to restore the agreement on the avoided crossing frequency 
while keeping the agreement on the large separation (model 
Sib). We therefore have M\ > M . 

4.5.2. Dependance of M on the helium abundance 

By modifying the value of the initial helium abundance in 
model So, we show that M linearly decreases with Yq (see 

Fig. nao. 

Like in the previous section, we try to explain this de- 
crease. We search for the new mass Mi that verifies condi- 
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Fig. 15. Evolutionary tracks in the HR diagram (left) and evolution of the central density p c (right) for models with the 
same mass but different values of a ov (a ov = 0, 0.08 and 0.15 from light gray to dark gray). The observed mean value of 
the large separation is represented by the dotted line. 



tion C when we increase the helium abundance to Y% >Yq. 
If we keep the same mass M (model S2a)> the mean molec- 
ular weight fi increases due to the increase of the helium 
content. It can be shown from homology relations that if we 
keep the mass constant and modify /i, the stellar luminosity 
scales as L oc /x 4 and the radius scales as R oc fi Z2 where 
the exponent z 2 depends o n the mode of energy genera - 
tion inside the star (see e.g. iKippenhahn k, Weigertl fl"990') . 
Subgiants are generally burning hydrogen through the CNO 
cycle and we then have Z2 — 0.6. In this case, combining 
these two relations with the fact that L oc R 2 T^ S , we ob- 
tain that T e s oc _L 5 7 . Increasing the helium abundance in a 
model while maintaining the mass constant induces a trans- 
lation of its evolutionary track in the HR diagram to the 
left along a line of slope about 5.7. Since the slope is bigger 
than that of the iso- Ai/ line, the models with higher helium 
abundance are closer to the TAMS when they reach Ai/ b s 
(green dashed line in Fig. IT3"t O . This effect is enhanced by 
the fact that when increasing Yoi the iso- Av line is trans- 
lated to the left. To reproduce the observed value of Across 
it is therefore necessary to decrease the mass and we have 
Mi < Mo (model Sab)- 
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Fig. 16. Variations of the optimal mass M with the over- 
shooting coefficient a ov . Two regimes can be distinguished: 
for a strong overshooting, the models are close to the TAMS 
and effect (i) dominates (a ov > 0.1) whereas for a milder 
amount of overshooting, effect (ii) takes over (a ov < 0.1). 



4.5.3. Dependance of M on the metallicity 

Fig. IT3"e shows that the mass M linearly increases with the 
abundance of heavy elements [Z/X]. 

To explain this increase, we consider a metallicity 
[Z/X] 3 > [Z/X]o and we search for the new mass M3. 
If we keep the mass Mo of model So (model S3 a ), the 
main effect of an increase of the metallicity is that the 
amount of bound-free absorption, which comes from met- 
als, is higher. This leads to an increase of the opacity. In 
order to study its effect, we have computed homology re- 
lations, following exactl y the s ame steps as in Chap. 20 of 
Kipp enhahn fc WeigertJ (fl990). but considering the opac- 
ity k as a basic parameter instead of the mean molecular 
weight [i. We found that if we keep the mass constant while 
increasing the opacity, the evolutionary track in the HR di- 



agram is translated to the right on a line of slope about 
3.6 (we note that the models indeed show a translation to 
the right, but with a somewhat larger slope a little below 
5, see purple dashed line in Fig.[TSf). Normally, this should 
bring the models closer to the TAMS when they reach the 
observed value of the large separation. However, we observe 
in Fig. [T3]that the \so-Av line is shifted to the right when 
increasing the opacity and the models with higher metal- 
licity are in fact further from the TAMS (purple full line in 
Fig.fTBT). We therefore need to increase the mass so that the 
frequency of the avoided crossing is correctly reproduced, 
i.e. M 3 > M (model S 3h ). 
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4.5.4. Dependance of M on the amount of overshooting 

We now focus on the influence of overshooting on the mass 
M, which is more complex than the influence of the stellar 
parameters previously studied. Starting from model So, we 
increased the amount of overshooting and determined the 
mass M. Fig. [TBI shows that M increases with a ov for low 
overshooting coefficients (a v < 0.06). Above this limit, 
the mass M starts to decrease as a function of a ov . This 
decrease eventually becomes linear for a ov > 0.1. We note 
that the shape of M(a ov ) does not qualitatively change if 
we consider initial models other than So (even though the 
value of the transitional a ov does vary). 

Core overshooting is not supposed to play a direct role 
since the models we here consider are in the PoMS stage 
and their convective core has vanished. However, it has an 
impact on the g-mode frequencies through its influence on 
the past evolution of the star and the chemical composition 
in the center. Adding overshooting has two opposing effects 
on the frequency of the avoided crossing: 

(i) First, by extending the mixed region associated to the 
convective core, overshooting increases the size of the 
hydrogen reservoir available for the nuclear reactions 
in the center. As a result, the star stays longer in the 
MS stage, as clearly appears in Fig. [T21 Consequently, 
for a given fixed mass, an increase of the amount of 
overshooting brings the star closer to the TAMS when it 
reaches the observed value of the large separation. Based 
on what has been said before, we expect the frequency 
of the avoided crossing to decrease when the amount of 
overshooting increases. 

(ii) But also, adding overshooting increases the size of the 
mixed core during the MS stage and therefore, when it 
reaches the TAMS, its helium core is bigger. The layers 
where hydrogen remains in the star are less deep and 
the temperature in these regions is thus smaller. The 
contraction of the star during the PoMS stage needs 
to be more severe to trigger nuclear burning in these 
shells. As a result, the central density increases faster 
with overshooting than it does without overshooting, as 
can be seen in Fig. 1151 Since the frequency of the avoided 
crossing is tightly connected to the central density (see 
Sect. l4.3Tll . we expect z/ cross to grow when a ov increases. 

In fact, these two effects play a different role, depending 
on the amount of overshooting. For a mild overshooting, 
the stars are evolved enough after the TAMS when they 
reach (Az/) b s so that the second effect prevails (see Fig. 
[15]). When increasing the amount of overshooting of model 
So from a ov = to a ov — 0.08, Across increases. This means 
that the stellar mass needs to be increased to ensure that 
condition C is satisfied. This accounts for the increase of M 
at low a ov . 

On the other hand, for higher amounts of overshoot- 
ing, the stars are very close to the TAMS when they reach 
(A^) bs- The central density has not had time to increase 
enough and the first effect is dominant. Fig. ll5l indeed shows 
that if we further increase the overshooting to a ov = 0.15, 
the central density decreases and so does the frequency of 
the avoided crossing. Therefore, the mass needs to be de- 
creased to satisfy condition C. 

We note that due to the competition between these 
two effects, certain masses lead to two different models 



Table 3. Ranges and steps adopted for the free parameters 
of the computed grids of models. 



Parameters 


Range 


Step 


Y 


0.24 to 0.28 


0.01 


t*conv 


0.48 to 0.72 


0.04 


[Z/X] (dex) 


0.04 to 0.14 


0.05 




0.00 to 0.20 


0.025 



fulfilling condition C: one with a mild overshooting away 
from the TAMS, and one with a larger a ov , close to the 
TAMS. For instance for the physics chosen in Fig. [TBJ a 
mass of M — 1.263 M@ can satisfy condition C both with- 
out overshooting (a ov = 0) and with an overshooting of 
a ov = 0.15. This observation shows that for subgiants in 
the early PoMS stage, knowing the frequency of an avoided 
crossing is not enough to efficiently constrain the amount 
of core overshooting in a star. 

5. Application to HD 49385: model optimization 

5.1. Computation of two grids of models 

We co mputed two grids of mod els, one assuming the mix- 
ture of iGrevesse k Noelsl (|1993H (fur hter noted GN93) and 
the other the more recent mixture of lAsplund et al.l {2005) 
(further noted AGS05). The grids are built with varying 
values of the mixing length parameter a conv , the helium 
abundance Yq, the metallicity [Z/X] and the amount of 
overshooting a ov . The ranges of studied values for each of 
these parameters as well as the chosen steps are reported in 
Table [31 For each point in the grids, an optimization such 
as the one described in Sect. 14.41 was performed to deter- 
mine the stellar mass M and age f which satisfy condition 
C. We then computed the x 2 function defined by Eq. [5] 

5.2. Results of the fit 

The first striking comment which can be made about the 
fit is that the representation of the total x 2 as a function of 
the amount of overshooting shown in Fig. [T7] presents two 
distinct minima, one for lower values of a ov (a ov < 0.05) 
and one for moderate values of a ov (~ 0.19). This leads us 
to identify four different families of solutions corresponding 
to the two different mixtures we considered and either low 
or high values of a ov . 

We then would like to derive a quantitative interpre- 
tation of the x 2 values we obtain for the models in each 
of these four families. This is possible since (i) the mea- 
sured dBta_caii_t)e_consMered_to have a Gaussian distribu- 
tion ([Appourchaux et ail 119981 showed that for simulated 
spectra with a resolution equivalent to four months of data, 
the distribution of the error bars for the central frequencies 
of the modes is close to a normal distribution) and (ii) for 
the four families of solutions, the models are roughly linear 
in their parameters in the regions defined by the uncer- 
tainties we will derive for the fitted parameters. This latter 
condition has been checked a posteriori. 

Under these conditions, the goodness of the fit can be 
estimated by the computation of the reduced x 2 defined as 

X?ed = TT^ (18) 
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Table 4. Values of the stellar parameters obtained from the computation of our two grids of models. 





GN93 


AGS05 


Qov 


low a ov 
n 0.5+ - 01 

u - uo -0.05 


high Q ov 
0.19 ±0.01 


low a ov 
0.00+ ' 01 


high a ov 
0.19 ±0.01 


M/Mq 
Age (Gyr) 

R/ Rq 

T e g (K) 
logff 

Qconv 


1.285 ±0.017 
4.90 ±0.13 

1.959 ±0.012 
5870 ± 40 

3.961 ± 0.002 
0.52 ±0.01 


1.231 ±0.012 
5.04 ±0.18 

1.932 ±0.011 
5790 ± 30 

3.946 ± 0.001 
0.52 ± 0.02 


1.264 ±0.013 
4.88 ±0.11 

1.947 ±0.007 
5940 ± 40 

3.960 ±0.002 
0.52 ±0.01 


1.210 ± 0.021 
5.10 ±0.18 

1.917 ±0.011 
6080 ± 60 

3.954 ± 0.002 
0.56 ± 0.03 


Minimum Xred 


2.48 


2.25 


2.11 


1.81 




Fig. 18. Left: Echelle diagrams of the best models with low overshooting (green) and with high overshooting (red). The 
squares represent i = modes, triangles £ = 1 modes and diamonds £ = 2 modes. The observations are represented by the 
black filled circles. Right: Evolutionary tracks in the HR diagram of the two models (green dashed line: low overshooting, 
red solid line: high overshooting). 
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Fig. 17. Values of the reduced % 2 for the models computed 
with the mixture of GN93 (blue empty circles) and with the 
mixture of AGS05 (red filled circles) as a function of the 
amount of overshooting. For clarity reasons, the symbols of 
models computed with GN93 were slightly shifted to the 
right. The best models for each considered value of a ov 
were linked by a dotted line (GN93) and by a dashed line 
(AGS05). 



Fig. 19. Values of the reduced x 2 as a function of the metal- 
licity [Z/X]. Symbols are the same as in Fig. [TT1 



where N is the number of measured data and P the num- 
ber of free parameters of the fit. More importantly, we can 
use the values of the x 2 function to determine a confidence 
interval for the fitted parameters. Suppose we want to de- 
termine the uncertainty in the amount of overshooting a ov 
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for instance. For each value of a ov in the grid, we determine 
the set of the other free parameters which minimizes the x 2 
func tion (the ob t ained \ 2 are linked by dotted lines in Fig. 

ress et all (|2002f) showed that the difference between 
these x 2 values and the minimum x 2 of the whole family of 
solutions Xmin defined as 

Ax 2 =X 2 ~ Xmin (19) 

is distributed as a x 2 with 2 degrees of freedom. This means 
that the " real" value of the parameter a ov has a probability 
of 68.3% (resp. 95.4%, 99.7%) to belong to the region where 
Ax 2 < 1 (resp. Ax 2 < 4, A% 2 < 9). The values of a ov for 
which Ax 2 = 1 thus give an estimate of the 1-a error in this 
parameter. This method is applied on all the parameters to 
determine the uncertainties which are reported in Table 0] 
for the four different families of solutions. 

We thus obtain a precise estimate of the fundamental 
properties of HD 49385. The stellar mass is constrained 
with a precision of about 1% for each of the four identified 
families. If we combine them all together, we find a mass 
of M = 1.25 ± 0.05 M Q for HD 49385 (corresponding to 
a precision of about 4%). Similarly, we obtain a radius of 
R = 1.94 ± 0.03 Rq for the star (precision of about 1.5%). 
The stellar age is also well constrained: from Table 21 we 
derive r = 5.02 ± 0.26 Myr (precision of about 5%). We 
obtain a log 5 of 3.954 ± 0.009 for HD 49385. We note that 
the results of the fit are consistent with the first estimates 
given in Table [TJ The level of precision we achieve here is 
to a great extent due to the detection of the £ = 1 avoided 
crossing in the spectrum of HD 49385 (see Sect. [5]). 

As expected, we find that we are also able to obtain 
information about the interior of HD 49385. As we men- 
tioned at the beginning of this section, we obtain strong 
constraints on the amount of overshooting in the star which 
can be either very low or moderate (a ov = 0.19 ± 0.01). In 
the latter case, the models are very close to the TAMS. We 
know that stars are expected to spend very little time in 
this evolutionary stage, which makes it quite improbable to 
observe them during this period. Although this point in it- 
self does not justify to ignore these solutions, it means that 
we should treat them with caution. The best models of the 
two families we obtained share very similar values of all the 
observables available for HD 49385, as shown in Fig. [THl so 
that we cannot discriminate between them. We note how- 
ever that the models show the existence of an £ = 2 mixed 
mode in the frequency range of the observations, which has 
different frequencies for the best low-overshooting models 
(between 973 and 984 /iHz) and for the high-overshooting 
ones (between 998 and 1002 /xHz). If such a mode could be 
detected in the observed spectrum, we would have a way 
of discriminating between the two scenarios. For complete- 
ness sake, we checked how the frequencies of the peaks 7t2 
and 7T3 of D10 compare with the £ = 3 frequencies in our 
models. In both cases (low and moderate a ov ), i"2 and 7T3 
appear between 2 and 3 /iHz higher than the theoretical 
values (i.e. between 3 and 5 a). 

An interesting feature of the fit is that the mixing length 
parameter a CO nv is well constrained for HD 49385 and cor- 
responds to a value which is significantly smaller than the 
solar one, regardless of the family of solution. Combining 
the families all together, we obtain a con v = 0.55 ± 0.04 to 
be compared to the solar value a con v = 0.64. To strengthen 
this result, we have checked that we reach similar con- 
clusions when using the traditional mixing length theory 



instead of the formalism of ICanuto et al.l (|l996f ) . In this 
case, the best agreement is reached for a mixing length of 
c*mlt ~ 1-4, significantly smaller than the solar calibrated 
value aMLT ~ 1-85. 

We also remark that the solutions are in favor of the 
abundances of heavy elements derived by AGS05. This can 
be clearly seen in Fig.[T5] This result is strengthened by the 
fact that even with the mixture of GN93, the best models 
are found when considering the lowest end of the range of 
the tested metallicities, i.e. [Z/X] = 0.04 dex. 

These results confirm that the seismic constraints avail- 
able for HD 49385 provide valuable information about the 
interior of the star. However, they raise a certain number 
of questions to which the grid-of-model approach in itself 
cannot answer. Why do we find two distinct families of so- 
lutions depending on the value of a ov ? How come do we 
manage to constrain the mixing length parameter a conv so 
well in this star? Why do smaller metallicities provide bet- 
ter fits than larger ones? Of what use are the information 
conveyed by the 1 = 1 ridge distortion on the coupling be- 
tween the cavities? We try to address these questions in the 
following section. 

6. An in-depth analysis of the results of the fit 

6.1. Sources of contribution to the x 2 function 

For the models computed in Sect. [5l we identified three 
main sources of contribution to the value of the x 2 function. 

The position in the HR diagram For each point of the grid, 
the mass and age are fixed to M and f, which determines a 
certain position in the HR diagram. With the abundances 
of AGS05, the values of T e g and L contribute only weakly to 
the x 2 of the best models (it represents only 0.3% of the 
total x 2 for our best model). The overall fit is thus 
fully consistent with the spectroscopic constraints. 
With the abundances of GN93, the fits point toward an 
effective temperature about 3-cr smaller than the measured 
one (see Tabled]) which results in a significant contribution 
to the x 2 value. 

The oscillation of the frequencies due to the T\ bump We 
have seen in Sect. 12.1.21 that a clear oscillation is observed 
in the large separation profiles. The models also show an os- 
cillation with a comparable period, which is caused by the 
bump of Ti in the helium second ionization zone (HIZ) at 
an acoustic depth of the order of thiz ~ 0.8. We note that 
another oscillation related to the base of the convection 
zone (BCZ) should also be present with a period around 
120 /iHz (tbcz ~ 0.5), however its amplitude is too low to 
be detected in the observations. The best agreement with 
the observed oscillation is obtained for thiz ~ 0.798, which 
is consistent with the estimate of the acoustic depth of the 
glitch we obtained in Sect. l2~i~2l (r glit ch = 0.76 ± 0.04). Fig. 
I2"0"k shows the comparison between the £ = large separa- 
tion profile of our best model and that of the observations 
(red curve). We can see that the agreement is good, except 
at high frequency where the period of the oscillation ap- 
pears to be smaller in the observations than in the models. 
However, the observations of p modes over a wider range 
of radial orders would be required to determine whether or 
not this decrease of the oscillation period is significant. 
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Fig. 21. Values of the quantity 6 (see text) for models obtained by varying one by one the stellar parameters of the 
best model without overshooting in Sect.[SJ the mixing length parameter (top left), the initial abundance of helium (top 
right), the metallicity (bottom left) and the amount of overshooting (bottom right). The filled circles indicate models 
computed with the abundances of AGS05 and the open squares with those of GN93. 



Even among the models which satisfactorily fit both the 
position of the star in the HR diagram and the oscillation 
of the mode frequencies due to the HIZ, some models yield 
high values of the \ 2 function. For instance, the model rep- 
resented in blue in Fig. (computed with a conv = 0.68, 
Y Q = 0.24, a ov = 0.10, [ZJX] = 0.14 dex and the mixture 
of GN93) reproduces the £ = large separations quite well 
but the I = 1 large separation profile shows some striking 
differences with the observations at low frequency, inducing 
a high value of the % 2 function (x 2 e d = 7.44) . This shows 
that the 1=1 ridge for this model is less distorted than it 
is in the observations. Based on the conclusions drawn in 
Sect. [31 this means that the coupling between the p-mode 
cavity and the g-mode cavity is too weak. 



The intensity of the 1=1 ridge distortion Since the inten- 
sity of the £ = 1 ridge distortion is a major contribution 
to the x 2 function (for instance for the model plotted in 
blue in Fig. [501 the three lowest frequency £ = 1 modes 
account for 80% of the \ 2 value) we tried to evaluate its 
dependence with the free parameters of our fit, in order to 
understand the results we obtained in Sect. [SJ For this, we 
first needed to quantify the intensity of the ridge distortion 



caused by the avoided crossing in our models. We chose 
to use the frequency difference 8 of the two modes which 
most closely surround the avoided crossing. To illustrate 
why these modes provide a good indicator of the ridge dis- 
tortion, their frequencies were circled in Fig. [7] for both of 
the studied cases (weaker coupling and stronger coupling). 
A low value of 8 indicates a weak coupling between the two 
cavities and a high value of 8, a stronger one. For HD 49385, 
we have 

8 = vi,i2 - v\,i\ (20) 

We insist on the fact that the value of 8, and in fact more 
generally the curvature of the £ = 1 ridge, can be used to 
estimate the strength of the coupling between the cavities 
only if the frequency of the avoided crossing in the models 
matches the one of the observations (i.e. v\ \\ = v 1Tl for 
HD 49385). This stresses the great interest of having access 
to the frequency of the mode which behaves mainly as a g 
mode in the avoided crossing (^1,11 in our case). It is also 
one of the reasons why we focused on finding all possible 
models satisfying this condition in Sect.|U 

We started from the best model obtained without over- 
shooting in Sect. [5] and we varied the stellar parameters 
one by one. The results are shown in Fig. 1211 Two striking 
remarks can be made about the obtained plots. 
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Fig. 22. Panel a: Propagation diagrams of models A (red solid lines) and B (blue dashed lines). The thick lines represent 
the Brunt- Vaisala frequency, the thin lines the £ = 1 Lamb frequency. The dotted line indicates the frequency of the 
avoided crossing (identical for both models). We then represent the central profiles of the nuclear reaction rate e to t(r) 
(panel b), of the gradient of mean molecular weight V/i (panel c) and of the temperature gradient V (panel d) for both 
models. The colored dotted lines indicate the boundaries of the evanescent zone for each model. 



First, we observe that the intensity of the I = 1 ridge 
distortion can account for most of the results of our fit. The 
quantity 6 decreases when the mixing length parameter in- 
creases, as can be seen in Fig. 12Tb . Only with low values 
of a C onv is it possible to reproduce the observed value of 
6 (6 ohs = 29.31 ± 0.33 fiRz, based on D10). This justifies 
the fact that our fit points toward a lower mixing length. 
Similarly, 6 decreases with increasing metallicity. We thus 
understand why the lower abundances of AGS05 (and more 
marginally the lowest metallicity with the abundances of 
GN93) provide a better fit to the observations. The ridge 
distortion appears to be relatively independent from Y$, 
and we found in Sect. [5] that Yq can hardly be constrained 
for HD 49385. Finally, the function 6 depends in a more 
complex way on the amount of overshooting. Two ranges 
of Q v satisfactorily reproduce 6 h s - very low amounts of 
overshooting (a ov < 0.05) and moderate amounts of over- 
shooting (a ov ~ 0.2). This corresponds exactly to what 
we obtained in Sect. [SJ All this establishes that the curva- 
ture of the I = 1 ridge due to the avoided crossing plays a 
crucial role in constraining the structure of the interior of 
HD 49385. 

The second remark which can be made is that the curves 
of 6(a conv ), S(Z/X) and 6(a ov ) in Fig. l2"Tlshow a remark- 
able anti-correlation with the curves of M(a conv ), M(Z/X) 
and M(a ov ) shown in Fig. ITBIandlTSl This is most striking 



in the case of the variations of 6 and M with the amount 
of overshooting. This clearly suggests that the coupling be- 
tween the cavities is inversely proportional to the stellar 
mass. This would explain why the mass of HD 49385 is so 
well constrained in this study. The next section is dedicated 
to propose a justification of this phenomenon based on the 
interior of the computed models. 

6.2. Why are the coupling between the cavities and the 
stellar mass anti-correlated? 

To understand the relation between the coupling and the 
stellar mass, we considered two of our models for which 
the £ = 1 ridge distortion is different: model A (computed 
with a conv = 0.48, Y = 0.24, [Z/X] = 0.09 dex using 
the abundances of AGS05) and model B (same as model 
A but with a con v = 0.72). Model A reproduces well the 
ridge distortion (6a — 29.51, 0.6-er agreement with the ob- 
servations) contrary to model B (6b = 27.60, 5.2-er agree- 
ment). Consistent with the results obtained in Sect. 14.5.11 
model B is more massive than model A (Mb = 1.37 M© 
and Ma = 1.20 Mq). We note that the results presented in 
what follows do not qualitatively change if we consider two 
other models with conflicting values of 6. 

Fig. I2"2"k shows the propagation diagrams for models A 
and B. The i = 1 Lamb frequencies of the two models al- 
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Fig. 20. Top: Profiles of the I = large separations for 
the observations (black filled circles and dotted lines) and 
for two models which reproduce well the oscillation due to 
the Ti bump. The red solid curve corresponds to the best 
model obtained in Sect. [5] and the blue dashed curve to a 
model computed with a conv = 0.68, Y Q = 0.24, a ov = 0.10, 
[Z/X] — 0.14 dex, using the GN93 mixture. Bottom: Same 
as the top panel for the 1=1 large separations. 



most overlap in the evanescent zone, but the Brunt- Vaisala 
frequency of model A in this region is larger than that of 
model B. This confirms that the coupling between the g- 
mode cavity and the p-mode cavity is stronger for model 
A. From Eq. [131 we know that the Brunt- Vaisala frequency 
depends on the temperature gradient V and on the gradient 
of mean molecular weight V/i. From the end of the main 
sequence phase to the present age, the profiles of these two 
quantities have been reshaped by the nuclear reactions in 
a shell. 

Since model B is more massive than model A, the tem- 
perature in its interior is larger, which causes the nuclear 
reactions of the CNO cycle to be more efficient in the hy- 
drogen burning shell. Fig. I2"2"b shows that the peak in the 
nuclear reaction rate £tot is indeed about twice as strong 
for model B. As a result, the hydrogen content in the layer 
where the reactions occur is more severely depleted and 
the corresponding peak in the gradient of chemical compo- 
sition is larger. It is therefore logical that in a small region 
above the hydrogen burning shell (which coincides with the 
evanescent zone), the gradient V/i, is smaller for model B. 
(see Fig. l22b). Based on Eq. [131 this suggests that the cou- 
pling between the cavities decreases when the stellar mass 
increases. Similarly, since nuclear reactions are more effi- 
cient in higher-mass stars, the ratio l/m (where l(r) is the 
luminosity going through the shell of radius r) is larger for 
model B and therefore the temperature gradient V is larger 
in the evanescent zone. This is confirmed by Fig. I22H . This 
effect also points toward a weaker coupling for higher-mass 
stars, according to Eq. [T31 
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Fig. 23. Propagation diagrams of two models without mi- 
croscopic diffusion (black dashed lines) and with diffusion 
(purple solid lines). The models both satisfy condition C 
and share the same values of the free parameters. The sym- 
bols are the same as in Fig. 



6.3. Remark on the effect of microscopic diffusion 

Microscopic diffusion was not included in the grids of mod- 
els presented in Sect. [SJ The main reason is that the time 
required to compute models including diffusion is several 
tens of times longer than when it is neglected, which would 
have severely limited the mesh of our grids. 

However, we have also computed a looser grid of models 
including microscopic diffusion in order to determine how 
this process influences our re sults. For this, we a dopte d 
the simplified formulation of iMichaud fc Proffittl ()l993h . 
Models were computed with either a low amount of over- 
shooting (a ov = 0) or a moderate one (a ov = 0.2), in ref- 
erence to the two families found in Sect. [5] We imposed an 
initial abundance of heavy elements such that the current 
surface metallicity matches the observed one. The mixing 
length parameter was varied from 0.48 to 0.72 (step of 0.04) 
and the initial helium abundance from 0.24 to 0.28 (step of 
0.02). Like before, for every set of parameters, we searched 
for a model satisfying condition C. 

We again obtained a best model offering a very satis- 
factory match to the observations (Xrcd = 2.1). All of its 
characteristics fall into the 1-<t error bars obtained with- 
out diffusion, except for the age which is smaller (4.1 Gyr 
for the best model) . The main reason why diffusion hardly 
changes the results is that it appears to have little influ- 
ence on the coupling between the cavities (which we have 
shown to play a dominant role in constraining the param- 
eters of the models). This can be seen in Fig. [23l which 
shows the propagation diagrams of two models computed 
with and without diffusion, both satisfying condition C and 
sharing the same values of the free parameters. We observe 
that the profiles of their Brunt- Vaisala frequency almost 
overlap. This can be understood as follows. The largest ef- 
fect of diffusion on the models is that the peak in the pro- 
file of V/i caused by the withdrawal of the convective core 
is smoothed. This should have an impact on the coupling 
since N oc V/i. However, this peak is rapidly reshaped by 
the nuclear reactions in shell. Since the center of the star is 
radiative, the nuclear reactions are quickly in equilibrium. 
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The abundances of the chemical elements in the core there- 
fore correspond to their equilibrium quantities which do not 
significantly vary when including diffusion. The nuclear re- 
action rates thus remain quite similar to the case without 
diffusion and so does the coupling between the cavities (see 
Sect.©. 

We note that interestingly, the solutions with high over- 
shooting are found less satisfactory when including micro- 
scopic diffusion. This probably has to do with the fact that 
the models with a ov — 0.2 which satisfy condition C are 
very close to the TAMS. Consequently, the gradient of the 
mean molecular weight has not yet been reshaped by the 
nuclear reactions in shell. A deeper analysis, out of the 
scope of the present paper, would be required to investi- 
gate whether or not the inclusion of microscopic diffusion 
can rule out the solutions with a moderate amount of over- 
shooting. 



7. Conclusion 

In this paper, we performed a detailed seismic modeling 
of the solar-like pulsating subgiant HD 49385. We first in- 
vestigated the possibility that the peculiarities observed in 
the oscillation spectrum of the star by D10 could be due 
to the existe nce of mixed modes. By exten ding a toy-model 
proposed bv lChristensen-Dalsgaardl (|2003l ). we showed that 
£ = 1 avoided crossings involve more than two modes and 
induce a characteristic and easily recognizable distortion of 
the £ = 1 ridge in the echelle diagram. We found PoMS 
(post main sequence) models with an I — 1 avoided cross- 
ing which can account for both the observed curvature of 
the ridge and the presence of one of the peaks detected 
outside the identified ridges. On the other hand MS (main 
sequence) models fail to reproduce these features. We thus 
established a firm detection of mixed modes in the spec- 
trum of HD 49385 as well as the PoMS status of the star. 
We also pointed out that the curvature of the I = 1 ridge 
depends on the coupling between the p-mode and g-mode 
cavities and should therefore bring information about the 
structure of the star in the evanescent zone which separates 
the cavities. 

We then remarked that the methods usually applied 
to model stars needed to be adapted to the special case 
of stars with avoided crossings. We therefore proposed a 
new approach to the grid-search modeling. We first showed 
that, for a given physics, the combined knowledge of the 
frequency of an avoided crossing and of the mean large 
separation is enough to obtain very precise estimates of 
the stellar mass and age (which we noted M and f in this 
paper). We then described a method to determine M and 
f in a systematic way in a grid of models. 

This method was applied to gain insights on how the 
stellar mass M varies with the different stellar parameters 
and efforts were made to physically understand these vari- 
ations. To model HD 49385, we used the proposed method 
to compute two grids of mod els, one assuming the h eavy 
element abundances of GN93 dGrevesse fc Noelslll993h an d 
the other assuming those of AGS05 ()Asplund et al.ll2005h . 
In this framework, we were able to strongly constrain the 
mass of HD 49385 (M = 1.25 ± 0.05 Mq) and its age 
(r = 5.02 ± 0.26 Myr). The stellar radius was found with a 
precision as good as 1% (R = 1.94 ± 0.03 Rq). 



We also obtained constraints on the physics of the in- 
terior of the star. Two different families of solutions were 
found, depending on the amount of overshooting that ex- 
isted above the convective core during the main sequence 
phase: one with a very low amount of overshooting (a ov < 
0.05), and the other with a moderate amount of overshoot- 
ing (a ov = 0.19 ± 0.01). The models of the latter fam- 
ily provide the best agreement with the observations, but 
they are very close to the TAMS and thus in a stage we 
should be unlikely to observe. The mixing length param- 
eter was found to be significantly smaller than the solar 
one (a CO nv = 0.5 5 ± 0.04 compared to a = 0.64 using 
the formalism of ICanuto et al.l|l996l ). Finally, we showed 
that the models comp uted with the revised heavy-element 
abundances of AGS05 (|Asplund et al.|[2005l ) provide a bet- 
ter ma tch with the observatio ns than the abundances of 
GN93 dGrevesse fc Noelsl[T99l . 

A more advanced study of the models of our grids 
enabled us to explain the results we obtained about the 
overshooting, the mixing length and the metallicity for 
HD 49385, in terms of stellar structure. We established 
that the intensity of the I — 1 ridge distortion associated 
with the observed avoided crossing plays a crucial role in 
constraining the parameters of the models computed for 
HD 49385. We showed that it can account for most of the 
obtained results. We also observed that the intensity of 
the coupling between the p-mode and g-mode cavities is 
strongly anti-correlated to the stellar mass in our models. 
We therefore suggested that these two quantities might be 
related to each other. By comparing the models of our grids, 
we managed to establish the link between them, mainly 
due to the fact that the temperature and chemical com- 
position profiles (upon which the Brunt- Vaisala frequency 
depends) strongly depend on the nuclear reaction rate in 
the hydrogen burning shell, and thus on the stellar mass. 
We therefore a posteriori understand why the mass is so 
well constrained in HD 49385. 

The subgiant HD 49385 is the first target for which a 
detailed modeling was led trying to reproduce all the prop- 
erties of an avoided crossing. It confirmed to a large degree 
that the detection of mixed modes provides unprecedented 
opportunities of probing the deep interior of stars. It also 
gave the opportunity to show that £ = 1 avoided cross- 
ing have a second interest: they enable us to probe the 
evanescent zone by yielding an estimate of the intensity of 
the coupling between the p-mode and g-mode cavities. We 
found here that this coupling strongly depends on the stel- 
lar mass and that we can thus obtain indirect constraints on 
the internal structure of the star. This shows that subgiants 
have the potential to contribute to the on-going debates 
which currently exist about the amount of overshooting at 
the boundary of convective zones, or the abundances of 
heavy elements in stars. This is all the more exciting since 
the space mission Kepler already claimed t he detection of 
£ = 1 avoided crossings i n several targets (jChaplin et al.l 
I20101 iMetcalfe et al.ll20i"oh . 
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Appendix A: Location of iso-Ai/ regions in the HR 
diagram 

We investigate the location in the HR diagram of mod- 
els which have the same large separation Av but different 
masses. Let us consider two such models with masses M 
and M'. Since the large separation is proportional to the 
square root of the mean density in the star, we have 




(A.l) 



It has been observationally established that there exists 
a very tight relation between the large separation and 
the frequency of m aximum power of the oscillations i/ maj( . 
IStello et all (|2009T) have found Av oc v^. We gather 
from this result that models which have the same large 
separation also sha re equivalent values of f ma x- Following 
I Brown et al.l ()1991l ). we expect v max to be proportional to 
the acoustic cut-off frequency and we have 

By combining Eg. IA.ll and lA.2[ we obtain R ~ T^ 2 . Since 
the luminosity is such that L ~ R 2 T^ S , we conclude that 
the models which have the same large separation satisfy 

L oc T c 5 ff (A.3) 

They are therefore located in the HR diagram on a line of 
slope equal to 5. This property is very well satisfied by the 
stellar models computed in this work. 



Gai, N., Basu, S., Chaplin, W. J., & Elsworth, Y. 2011, ApJ, 730, 63 
Gelly, B., Lazrek, M., Grec, G., et al. 2002, A&A, 394, 285 
Grevesse, N. & Noels, A. 1993, in Origin and Evolution of the 

Elements, ed. N. Prantzos, E. Vangioni-Flam, & M. Casse 

(Cambridge University Press, New York), 15—25 
Kippenhahn, R. & Weigcrt, A. 1990, Stellar Structure and Evolution, 

ed. (Springer, Berlin) 
Kjeldsen, H., Bedding, T. R., & Christensen-Dalsgaard, J. 2008, ApJ, 

683, L175 

Kjeldsen, H., Bedding, T. R., Viskum, M., & Frandsen, S. 1995, AJ, 
109, 1313 

Koch, D. G., Borucki, W. J., Basri, G., et al. 2010, ApJ, 713, L79 
Lebreton, Y., Monteiro, M. J. P. F. G., Montalban, J., et al. 2008, 
Ap&SS, 316, 1 

Metcalfe, T. S., Monteiro, M. J. P. F. G., Thompson, M. J., et al. 
2010, ApJ, 723, 1583 

Michaud, G. & Proffitt, C. R. 1993, in IAU Colloq. 137: Inside the 
Stars, ed. W. W. Weiss & A. Baglin (Astronomical Society of the 
Pacific Conference Series, San Francisco), Vol. 40, 246—259 

Miglio, A. & Montalban, J. 2005, A&A, 441, 615 

Morel, P. 1997, A&AS, 124, 597 

Morel, P., van't Veer, C, Provost, J., et al. 1994, A&A, 286, 91 
Mosser, B. & Appourchaux, T. 2009, A&A, 508, 877 
Osaki, J. 1975, PASJ, 27, 237 

Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 

2002, Numerical recipes in CH — h : the art of scientific computing, 

ed. Cambridge University Press (New York) 
Provost, J., Berthomieu, G., Martic, M., & Morel, P. 2006, A&A, 460, 

759 

Roxburgh, I. W. 2009, A&A, 506, 435 

Samadi, R., Kupka, F., Goupil, M. J., Lebreton, Y., & van't Veer- 

Menneret, C. 2006, A&A, 445, 233 
Scuflaire, R. 1974, A&A, 36, 107 

Scuflaire, R., Montalban, J., Theado, S., et al. 2008, Ap&SS, 316, 149 
Stello, D., Chaplin, W. J., Basu, S., Elsworth, Y., & Bedding, T. R. 

2009, MNRAS, 400, L80 
Tassoul, M. 1980, ApJS, 43, 469 

von Neuman, J. & Wigner, E. 1929, Zhurnal Physik, 30, 467 



References 

Aerts, C, Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, 

Asteroseismology, ed. Springer (Berlin) 
Aizenman, M., Smeyers, P., & Weigert, A. 1977, A&A, 58, 41 
Angulo, C, Arnould, M., Rayet, M., et al. 1999, Nuclear Physics A, 

656, 3 

Appourchaux, T., Gizon, L., & Rabello-Soares, M. 1998, A&AS, 132, 
107 

Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in Cosmic 
Abundances as Records of Stellar Evolution and Nucleosynthesis, 
ed. T. G. Barnes III & F. N. Bash, (Astronomical Society of the 
Pacific Conference Series, San Francisco), Vol. 336, 25 

Baglin, A., Auvergne, M., Barge, P., et al. 2006, in The CoRoT 
Mission, ed. M. Fridlund, A. Baglin, J. Lochard, L. Conroy (ESA 
Special Publication, Noordwijk), Vol. 1306, 33 

Barban, C, Michel, E., Martic, M., et al. 1999, A&A, 350, 617 

Brown, T. M., Gilliland, R. L., Noyes, R. W., & Ramsey, L. W. 1991, 
ApJ, 368, 599 

Canuto, V. M., Goldman, I., & Mazzitelli, I. 1996, ApJ, 473, 550 
Carrier, F., Eggenberger, P., & Bouchy, F. 2005, A&A, 434, 1085 
Chaplin, W. J., Appourchaux, T., Elsworth, Y., et al. 2010, ApJ, 713, 
L169 

Christensen-Dalsgaard, J. 1981, MNRAS, 194, 229 
Christensen-Dalsgaard, J. 2003, Lecture Notes on Stellar Oscillations, 
5th edn. 

Christensen-Dalsgaard, J. & Thompson, M. J. 1997, MNRAS, 284, 
527 

Deheuvels, S., Bruntt, H., Michel, E., et al. 2010, (D10), A&A, 515, 
A87 

Deheuvels, S. & Michel, E. 2010, Ap&SS, 328, 259 

Di Mauro, M. P., Christensen-Dalsgaard, J., Kjeldsen, H., Bedding, 

T. R., & Paterno, L. 2003, A&A, 404, 341 
Di Mauro, M. P., Christensen-Dalsgaard, J., Paterno, L., & D'Antona, 

F. 2004, Sol. Phys., 220, 185 
Dziembowski, W. A. & Pamyatnykh, A. A. 1991, A&A, 248, Lll 
Gabriel, M. 1980, A&A, 82, 8 



